SFB/CPP-10-111 
TTP 10-47 



Towards Higgs boson production in gluon fusion to 



We consider the Higgs boson production in the gluon-fusion channel to next- 
to-next-to-leading order within the Minimal Super symmetric Standard Model. In 
particular, we present analytical results for the matching coefficient of the effective 
theory and study its influence on the total production cross section in the limit 
where the masses of all MSSM particles coincide. For supersymmetric masses below 
500 GeV it is possible to find parameters leading to a significant enhancement of 
the Standard Model cross section, the X-factors, however, change only marginally. 
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1 Introduction 



The hunt for the Higgs boson is performed with great enthusiasm both at the Fermilab 
Tevatron and at the CERN Large Hadron CoUider (LHC). The Higgs boson constitutes 
the last missing piece in the Standard Model (SM) of particle physics and could explain 
the generation of mass terms in the Lagrange density. The latter is also true in many 
extensions of the SM, such as supersymmetry. 

In the recent years many higher order corrections to production and decay of the Higgs 
boson have been computed. The leading order (LO) and the next-to-leading order (NLO) 
corrections to the production cross section in gluon fusion have been considered in Refs. 
[iHl]. The most recent computations within the SM that appeared in the end of 2009 [SH^ 
describe the Higgs boson production in gluon fusion at the next-to-next-to-leading order 
(NNLO) taking the finite top quark mass into account. An important outcome of these 
papers is the justification to use the effective theory (built for the infinite top quark mass) 
for the evaluation of the cross section as it was done in Refs. P4T0]. As long as one re- 
stricts the analysis to the top quark/squark sector it is suggestive that this conclusion also 
holds for the Minimal Supersymmetric Standard Model (MSSM) since the correspond- 
ing calculation would have a similar structure. At NLO this has been demonstrated in 
Refs. [II1[I2]- In this paper we assume that it also holds at NNLO in situation when all 
supersymmetric particles are heavier than the top quark. Note that we explicitly exclude 
the bottom quark/squark sector from our consideration since in that case the effective- 
theory approach as presented in this paper is not apphcableQ Thus, our results only 
apply to small and moderate values of tan/3, the ratio of the vacuum expectation values 
of the two Higgs doublets, where the contribution of the bottom quark/squark sector is 
numerically negligible. 

The supersymmetric NLO correction to the gluon fusion process has been considered 
in various papers. In Refs. [HI [15] the Higgs boson production cross section and the 
decay rate have been computed in the effective theory framework and thus apply to Higgs 
boson masses less than about 200 GeV. The results have been confirmed by a (numerical) 
calculation in the full theory [12] (see also Refs. [TB|[T7]). It should be mentioned that 
the squark contribution has been considered separately in Refs. [TSKin]. Furthermore, an 
estimate of the NNLO supersymmetric corrections was presented in Ref. [20] (see also 
Ref. [11]) assuming that the NNLO corrections to the matching coefficient in the MSSM 
and the SM are numerically of the same order. 

In this paper we consider the strong interaction part of the MSSM and evaluate the total 
production cross section of an intermediate mass Higgs boson in the process of gluon 
fusion. As motivated above, we employ the effective theory obtained in the formal limit 
where all supersymmetric particles and the top quark are infinitely heavy, which leads to 

^See, e.g., Refs. [13l[T4] where the bottom quark/squark sector is discussed at NLO in detail. 
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the Lagrange density 



with 



Ol = \g%G''^\ (2) 

Ci is the coefficient function containing the remnant contributions of the heavy particles 
and is the effective operator. jCqIj^ is the QCD Lagrange density with five active 
flavours. In our case (p^ denotes the light CP-even Higgs boson of the MSSM. is the 
gluonic field strength tensor of the standard five-flavour QCD and the superscript "0" 
indicates bare quantities. Note that the renormalization of 0°/ is of higher order in the 
electromagnetic coupling constant. 

The renormalization of Oi is discussed in detail in Refs. [211 122]. Fo^' completeness we 
reproduce here the renormalization constant which in the MS scheme is given by 

0, = ZuOl Zn=(l--—] , (3) 

where ci = 4 — 2e is the space-time dimension and Z?^^-* = fi^^^[oL^ is the /3 function of the 
standard five-flavour QCD with /J^^) = -{af jT^f^f , , . and = 11/4 -n//6. Since 
C\0\ = CiOi, the coefficient function is renormalized with l/Zu. 

Following Ref. [23], we slightly redefine the coefficient function and the operator in order 
to have objects which are separately renormalization group invariant. Introducing 

= ^C,, (4) 



with 



we may choose independently the "soft" and the "hard" renormalization scales. Our 
following analysis of the Higgs production cross section takes advantage of those properties 
of Og and Cg. Note that the LO piece of Cg is not proportional to ag- 
in this paper we consider the strongly interacting part of the MSSM and assume for 
simplicity the degenerate supersymmetric mass spectrum. The evaluation of the coeffi- 
cient functions is performed within Dimensional Reduction (DRED) [23] where we follow 
Refs. [2nil2S] for the practical implementation. 
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The remainder of the paper is organized as follows: in the following section we briefly 
revisit the SM and in particular discuss the effect of B^^^ introduced above. We further- 
more evaluate the SM coefficient function within DRED which is a useful preparation for 
the SUSY QCD calculation. Next, in Section [3] we discuss the three-loop calculation of 
the MSSM coefficient function which is used in Section H] to compute the cross section to 
NNLO. Section |5] contains the summary and conclusions. 



2 Revisiting the SM 

As mentioned in the Introduction, the gluon fusion process has been discussed extensively 
in the literature. The purpose of this section is twofold: on the one hand we want to discuss 
the effect of B^^^ introduced in Eq. (jl]) on the numerical predictions, on the other hand we 
repeat the SM calculation in the framework of DRED as a preparation for the calculation 
in the supersymmetric theory. 



2.1 Separation of hard and soft scales 



The SM coefficient function expressed in terms of the on-shell top quark mass and the 
five-flavour strong coupling read^ [22l[27] (see also Ref. [28]): 



SM 



af^ 11 




19 

16 M? 



+ ni 



67 1 , /i^ 
h - In 

96 ^ 3 M? 



(6) 



where ri; = 5 is the number of massless quarks. 

In Fig. [m^a) we show the production cross section a{pp H + X) for the LHC with 
the center-of-mass energy of 14 TeV at LO, NLO and NNLO for Mt = 173.3 GeV 
computed with MSTW08 [30] parton distribution function (PDF) set. 

The chosen PDFs fix the values of af\Mz) to LO, NLO and NNLO accuracy used as the 
starting point in order to obtain the strong coupling for the other choices of the renormal- 
ization scale. The bands are obtained by varying the renormalization and factorization 
scales fJLr = l^f = 1^ between Mh/^ and Mh [311l32j^ with n = Mh/2 as central value. For 
comparison we present the NNLO result computed with ABKM09 PDFs [M] (the dashed 
line). It is remarkable that the difference in the central values for the different PDFs is 
comparable to the uncertainty due to the scale variation. 



^Here and in the following we suppress the ^-dependence of a^. If not stated otherwise we have 
as = as (fi) ■ The same is true for the DR mass parameters rrit and msusy introduced in Section |3l 
•^For the other opinions concerning the scale variation see Ref. [55] , 
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Figure 1: (a) The SM cross section as a function of Mh- The LO, NLO, and NNLO results 
are represented by the lower, middle (coarse filling) and the upper bands, respectively. 
Band width is due to simultaneous variation of renormalization and factorization scales. 
The dashed line corresponds to the NNLO prediction using ABKM09. (b) Ratio R = 
a{fis, Mt) /a{fis, fig) for = 120 GeV as a function of fig- Dotted, dashed (barely 
visible) and solid lines correspond to LO, NLO and NNLO, respectively. 
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The three-loop term of Eq. contains the logarithm \n{fi^/M^) which can potentially 
lead to numerically big effects if the scale /i assumes very small or very large values. With 
the help of Cg and Og introduced in Eq. (j4]) we may introduce separate scales fih and 
fig for the hard and the soft parts of the cross section a = a^fisyfJ^h)- The dangerous 
logarithms can then be avoided by choosing fih ~ Mt which is the natural scale for the 
coefficient function, and varying fig around as done above. Note that for /i^ = fi^ 

the quantity i?^^^ drops out of the analysis. 

In Fig. [It^b) we plot the ratio R = a{fis, Mt) /a{fis, ^g) for = 120 GeV. In the nu- 
merator the hard scale is kept fixed while in the denominator we use the common ap- 
proach P4T0] and identify the scales as in Fig. [U^a). At LO i? = 1 by construction. The 
NLO corrections (dashed curve which is hardly visible) are quite small and amount to 
less than ±2% which is only slightly modified after inclusion of the NNLO results (solid 
line). From these considerations we conclude that the resummation of ln(yU^/M^^) terms 
has only a small numerical effect. On the other hand, the framework of Eq. (jl]) is quite 
useful in the context of supersymmetric corrections discussed below since it allows to fix 
the hard scale fit and vary /i^. 

2.2 Cf^ in DRED 

Whereas DRED is the preferable framework for the two- and three-loop calculations in a 
supersymmetric theory, in the SM the application of DRED introduces several so-called 
evanescent couplings which can make practical calculations quite tedious. This has been 
discussed extensively in the literature (see, e.g., Refs. [35l - [39] ). Nevertheless, for our 
purpose it is useful to evaluate the QCD corrections to the matching coefficient using 
DRED since there are several new features which also apply to the MSSM. 

We implement DRED with the help of the so-called e scalars. They have mass dimension 
(4 — d) and complement the ci- dimensional gluon field so that the gauge fields and the 
fermions have exactly four degrees of freedom. In a practical calculation the e scalars 
induce new Feynman rules (see, e.g., Ref. [ID]) which (in non-supersymmetric theories) 
involve new couplings that have to be renormalized independently from the gauge cou- 
pling. 

In our matching calculation it is convenient to define Cf^ in such a way that the effective 
theory is built from the five-fiavour QCD fields regularized with DREG. As a consequence, 
we have to integrate the e scalar field out together with the top quark in order to guar- 
antee the transition from DRED to DREG in the matching procedure. This means that 
formally ^ and even ^ Mu since the latter is identical to zero for the matching 
calculation. In the end, we take the limit — )■ in agreement with our renormalization 
condition (see below). 

We evaluate the matching coefficient in the limit Mt ^ M^ ^ and apply the rules of 
asymptotic expansion together with a naive expansion in the external momenta in order 
to project out the Lorentz structure of the Higgs-gluon-gluon vertex (see the following 



6 



section for details). Consequently, there are (sub-) diagrams which contain only e scalars 
as massive particles. Due to the dimension of the underlying integrals they lead to terms 
proportional to and in the sum of the bare diagrams. Thus it is necessary to 
keep positive powers in in counterterms entering the calculation. 

The top quark mass and the strong coupling need to be renormalized at two-loop order. 
Furthermore, one has to make sure that the external gluon fields of the full and the 
effective theory Green functions are identical. There are two operators involved in the 
renormalization procedure which deserve additional attention: the e-scalar mass and the 
coupling of the Higgs boson to e scalars. The corresponding part of the Lagrange density 
is 



1 



(7) 



As mentioned above, we keep ^ for the matching procedure, however, in the final 
expression for Cf^ we require that the on-shell vanishes. This leads to the well- 
known counterterm (see, e.g., Refs. [2S1I1I]) which we need to one-loop accuracy. For our 
purposes we have to keep the terms proportional to and restore the dependence on 
all evanescent couplings. 

The second term in Eq. ([7]) is not present at tree-level, however, quantum corrections 
induce such a structure [12]. We require that the renormalized parameter vanishes 
since we want to avoid a non-zero coupling between the Higgs boson and the unphysical e 
scalars. As a consequence of this condition, we have to introduce a non-zero bare coupling 
of the Higgs boson to e scalars. Furthermore, in the counterterm contribution to A^ we 
only have to keep terms not proportional to A^. We refrain from listing the explicit results 
which become quite lengthy due to the evanescent couplings. 

It is an important check of our calculation that the negative powers in cancel after 
the renormalization of A^ and M^. 

Our final result for the coefficient function reads: 
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where It = ln[/i^/(mf(yU^))^ 



a 



DR 



is the strong coupling with six active flavours in the 



DR scheme and is the evanescent coupling originating from the e scalar-quark vertex. 
Note that the evanescent coupling stemming from the four-e- vertex cancels in Eq. ([8]). 
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The transition to the five-flavour QCD is achieved with the help of 



„(5) r 



.-1 



+ 



(5)' 

as 

TT 
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11 13, 1 ,2 

^ H 

9 24 36 * 
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h fii — 

36 12 



(9) 



From Eqs. (|H]) and we recover Cf^"^ of Eq. (iniu which is a welcome cross check on 
all individual steps of our calculation, in particular on the counterterms originating from 
Eq. ([7]). In the next section the same framework is applied to the MSSM. 



3 Coefficient function in the MSSM to three loops 

In this section we describe the evaluation of the coefficient function Ci of Eq. ([1]) within 
the strong sector of the MSSM. Note that the operator Oi contains the effective inter- 
action of two, three or four gluons with the Higgs boson and we are free to choose any 
corresponding Green functions in order to extract Ci. We decided to consider the two- 
gluon-Higgs boson amplitude since it involves fewer diagrams. We need to keep the two 
gluon momenta different until the Lorentz structure of Oi responsible for the two-gluon 
coupling is projected out. 

At one-, two- and three-loop order 10, 671 and 49632 Feynman diagrams have to be 
considered. Some typical diagrams are shown in Fig. [2l In addition to the particles of 
supersymmetric QCD we also have the e scalar. Higgs boson couplings to the light quarks 
and the corresponding squarks are omitted. Furthermore, we neglect the contribution of 
the Higgs boson-squark coupling suppressed by M|/m^ (for the Feynman rules see, e.g.. 
Appendix A of Ref. [TT]l 

Our calculation is based on a chain of programs designed to minimize the manual input. 
First, we generate all the Feynman diagrams with QGRAF [32], then transform the output 
to FORM [13] notation with the help of q2e and apply asymptotic expansion (see, e.g., 
Ref. [B]) using exp [13116]. 

Due to the gauge invariance the Higgs-gluon-gluon coupling has the form 

^,u,ab ^ A6'''{g^''q,-q2-q'[q^) , (10) 

where A is a scalar function depending on = 2gi ■g2, M^, /i^, and the mass parameters 
of SUSY QCD. qi and q2 are the outgoing momenta of the gluons with polarization vectors 
e^{qi) and e'^{q2) and colour indices a and b. We begin by convoluting the vertex Feynman 

^Note that the functional form of Eq. ^ is independent from the quark mass definition. 
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Figure 2: Sample diagrams contributing to Ci at one, two and three loops. The symbols 
t, ti, g, g, and e denote top quarks, top squarks, gluons, gluinos, Higgs bosons and e 
scalars, respectively. 

diagrams with the projector 



^^^^^^'^^^ - ^ 8id-2)iq,.q,r ' ^^^^ 



SO that only scalar integrals remain. After the Taylor expansion to the linear order in 
qi and q2 the factor (gi ■ ^2)^ in the denominator of Eq. (11 II) cancels and both momenta 
can be set to zero. Finally, we evaluate vacuum integrals which (after the asymptotic 
expansion) only contain a single scale. This is done using the package MAT AD 

Due to the transverse structure of the Higgs-gluon-gluon amplitude ^/^''•"^ an alternative 
projector to Eq. ffTTl) can be applied to obtain Ci. 



^'^'^^^^'^'^ " ^ %{d-2){q,.q,f • ^^^^ 



We have checked that in both cases we obtain the same result. 

One may consider different hierarchies between the mass scales involved in this calculation. 
As the first choice, we identify all supersymmetric masses with msusy and thus have two 



9 



expansion parameters, Xts = mt/rrisusy and Xst = M^/rrit. At three-loop order we have 
computed terms through O (x^g) and As a cross-check, we performed a similar 

calculation where the top quark mass and msusy were identified from the very beginning. 
Comparing the latter result to the former expansion in xts evaluated at Xts = 1 we find 
agreement to within a few per cent, which serves as a convergence test of our results. 

We choose to renormalize all parameters in the DR scheme In particular, we need 
the counterterms for a^, the top quark, and the top squark masses to two- loop, and 
the gluino and other squark mass counterterms to one-loop order. Furthermore, the 
counterterm for the mixing angle 6t between the top squarks is also required to the 
two-loop order. Although 6t is zero in our approximation, its renormalization in the non- 
degenerate one- and two-loop expressions leads to non-zero contributions at three-loop 
level. The corresponding counterterm can either be obtained from the exact one- and 
two- loop result for Ci [11], or (as done in this paper) via expansion in the mass difference 
between the two top squarks. Note that the same effect is also present in the three-loop 
calculation of the Higgs boson mass within the MSSM, as discussed in Ref. [26lll8]. All 
necessary counterterms can be found in the Appendix of Ref. [2H] except for the decoupling 
constant relating the gluon field in the SUSY-QCD to that in the five-fiavour QCD, given 
by 



AO 
^3 



(SQCD) 



1, 1 - 

+ fT'l—l'susy 

6 12 




n, —I, 



susy 



(13) 



where ^usy = ln[/i^/(rrasusy(/^^))^] and ai^^'"^^ denotes the DR coupling where all particles 
of the theory contribute to the running. In Eq. (fT3|) ri; = 5 counts the number of massless 
quarks and at the same time the number of squark fiavours different from top squarks. 
Similarly to the SM case in Section [221 we keep 7^ in the intermediate steps in order 
to match to the effective theory of Eq. ([T]). Thus, and the coupling A^ (cf. Eq. ([7])) 
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have to be renormalized with the help of the counterterms 
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(14) 



with U = ln[fi'^/M^], and 
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(15) 



where Cq, = cos a, S/? = sin/3, to, = tana, t/3 = tan/3 and x^t = fisusy/^t- Here a is 
the mixing angle between the weak and the mass eigenstates of the neutral scalar Higgs 
bosons, tan /3 is the ratio of the vacuum expectation values of the two Higgs doublets, 
and //susy is the Higgs- Higgsino bilinear coupling from the super potential. Note that 
the counterterm for and the three-loop vertex diagrams have been evaluated assuming 
the same hierarchy of the top quark and the top squark masses. Let us mention that 
the counterterm in Eq. fll5p is finite at one-loop order and only contains a single pole at 
two loops whereas the corresponding expression in QCD has a single and a double pole, 
respectively. 

In the first step we obtain Ci in terms of It is, however, more convenient to 

express the final result in terms of the five-flavour coupling in the MS scheme a 



(5) 



a 



(5),MS 



. The relation establishing this transition is as follows: 



a 



(5) 



(SQCD) 



(16) 



Two-loop corrections to the decoupling constant within supersymmetric QCD have 
been computed in Refs. [ISHSI]- Since (as i^i Eq. (ITB]) combines the decoupling of the 



heavy particles and the transition from DR to MS, we present the explicit result in the 
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limit of equal supersymmetric masses: 

JSQCD) / ^ ^ ^ 
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(17) 



All occurring mass parameters and the top squark mixing angle have been renormalized 
in the DR scheme. 



Finally, we present the coefficient function Ci expressed in terms of DR masses in the full 
theory, expanded in af^: 
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The analytic expression including 0{x^g) and x^^) is available on request from the 

authors. It is noteworthy that in the limit mgusy — )■ C)0 we recover exactly the SM result 
of Eq. ©. 



4 Higgs boson production at NNLO 

In this section we discuss the numerical effect of the corrections presented above. The 
effective-theory parts in the MSSM and the SM are identical, and thus we can borrow the 
framework of Ref. [6]. In addition, we incorporate the factor 5^^^ (in the MS scheme) as 
discussed in the Introduction. 
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On the full-theory side we consider the renormalization group invariant combination Cg = 
Ci/B^^'^ where both Ci and B^^^ are parameterized with the five-flavour strong coupling 
in the MS scheme. For the numerical evaluation of ^ we fix the (hard) renormalization 
scale at the on-shell top quark mass value, i.e. /x/i = Mf, which means that the DR mass 
parameters in Eq. (fTSl) have to be interpreted as msusy(^i) and mt{Mt). Any variation 
of /i/i leads to a non-degenerate supersymmetric spectrum and thus renders Eq. (|T8i) not 



applicable. As shown in Section |2] (Fig. E^b)), the variation of Hh only leads to a very 
small effect and thus we keep it fixed. 

Following the common practice, we write the hadronic cross section for the individual 
channels a^ 



(yhk[s] 



(19) 



where hk G {gg, qg, qq, qq, qq'} (g and q' denote different light quark flavours) and ao 
contains the exact dependence on the heavy masses at the LO. In the SM this is only the 
top quark mass, while in the MSSM also the top squark masses m^^ = m 
be taken into accountI§: 



''t2 



TTlsusy must 



Cq, 

2 /T\j2 
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(20) 



where n = 'imf/Mfj, Ts = Aml^^^/Mjj and 



m 



arcsin^(l/v^). 



t > 1 , 

t < 0. 



(21) 



The quantity T^hkis) in Eq. ( !T9|) . expanded in ah\^s)-, contains the convolutions of the 
PDFs with the n-loop {n = 0, 1, 2) effective-theory partonic cross sections. In the SM, 
such effective theory approach provides an excellent approximation at NNLO [5H7]. Once 
the LO cross section which depends on af'\fis) is factored out, we are left with a product 
of the residual hard coefficient depending on ai^\nh) and the remaining part of T^hkis) 
which is parameterized in terms of af\fis). We expand these factors and consistently 
discard higher orders in the strong coupling, decomposing the hadronic cross section into 
the LO, NLO and NNLO pieces. Introducing 

Cg = l + ^^c^^^+(^^\cf^ + ..., 




■'The subscripts hk mark the partons in the initial state. Below we sum over all possibilities. 
^The terms of order M^/rn^ are not displayed but are accounted for in the numerical evaluation. 
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Figure 3: LO (lower), NLO (middle) and NNLO (upper curves) Higgs production cross 
section as a function of mgusy The solid and dashed lines correspond to the MSSM and 
SM, respectively. 



we write the production cross section as 
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(23) 



where superscripts on Cg or S denote the loop order of the corresponding quantity. 
As input for the numerical evaluation we use the on-shell top quark mass Mt 



173.3 GeV [2H] and convert it to the DR scheme using two-loop SUSY-QCD correc- 
tions [52]. By default we use LO, NLO and NNLO PDFs by MSTW2008 collabora- 
tion [30]. The choice of the PDFs determines the values of a^^\Mz) at the considered 
order. Using the appropriate MS beta function, we then derive af'\jj,s) and af\ijLh). 

In Fig. |3]we show the total cross section as a function of a common supersymmetric mass 



scale rrisusy up to NNLO with tan/3 = 5, ii., 



susy 



200 GeV, vanishing mixing angle for the 
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top squarks and a pseudoscalar Higgs boson mass Ma = 1 TeV. The Higgs boson mass 
is computed with the help of the program H3m [26] as a function of mgusy to three-loop 
accuracy. In the range of mgusy varying between 250 GeV and 2 TeV Mh assumes values 
between 92 GeV and 117 GeV. Although experimentally excluded if the SM limits on Mh 
are assumed, we show the curves for illustration in the mentioned region for mgusy The 
mixing angle a is relatively stable so that ~ 0.99. The SM results are shown with 
the dashed lines. They develop mgusy dependence since we identify the SM Higgs boson 
mass with the MSSM value. Note that for msusy ^ 1500 GeV the corresponding Higgs 
boson masses are excluded by LEP [S3]. Nevertheless we continue the dashed curves to 
'^susy = 250 GeV to have a comparison with the MSSM curves. 

One observes that for the large values of msusy the MSSM curves approach the SM result 
whereas for mgusy ^ 400 — 500 GeV a significant enhancement of the MSSM cross section 
is observed. 

In order to determine the effect of higher order corrections it is convenient to define the 
so-called i^-factor as the ratio of the higher order cross section and the Born cross section 
with the soft scale fixed at /i^ = (M///2)^, 

The values of the i^-factor corresponding to Fig. [3] vary only slightly with msusy and are 
^NLO ^ ]^_5^ ^NNLO ^ pjQgg ^Yieii SM counterparts. 

We refrain from discussing the renormalization scale dependence since it is essentially 
inherited from the SM corrections discussed in Section [2l 



Our following scenario is based on the so-called "no-mixing" scheme as defined in Ref. 
Keeping tan/3 = 5 fixed, we vary Ma, the mass of the pseudoscalar Higgs boson, between 
200 GeV and 1 TeV, and compute the scalar Higgs mass and the supersymmetric particle 
masses with SOFTSUSY [55]. Averaging the masses of squark, stop and gluino, we obtain 
msusy ~ 1970 GeV. In Fig. H] we present the cross section computed in this scenario. Here 
the Higgs boson mass varies from 114 GeV to 117 GeV, and that affects the SM cross 
section (the dashed lines). Compared to the previous plot, the variation of a is more 
pronounced, reaching ?a 0.95 at the left boundary, while for the larger values of Ma 
this ratio approaches 1. The genuine supersymmetric corrections (arising from terms of 
(!?(l/msusy) in Eq. f llSp ) are small. Again, the /^-factors are very stable in the whole range 
of Ma and are K^^^ = 1.6 and i^NNLO ^ ^ g_ 



5 Conclusions 

In this letter we have computed supersymmetric NNLO SUSY-QCD corrections arising 
from the top quark/squark sector to the matching coefficient determining the Higgs-gluon 
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Figure 4: LO (lower), NLO (middle) and NNLO (upper curves) Higgs production cross 
section for the "no-mixing" scenario as a function of Ma- The solid and dashed lines 
correspond to the MSSM and SM, respectively. 



coupling in the effective theory. Regularization was done in DRED with the help of the 
so-called e scalars which made the matching to the five-flavour QCD very delicate. In 
particular, we had to renormalize the mass of the e scalar and its coupling to the Higgs 
boson at one- and two-loop orders, respectively. 

The approximation investigated in this paper, assuming zero mixing in the stop sector, 
produces significant effects only when mgusy is small, -C 1 TeV. However, the effect of 
zero mixing on Higgs mass together with the current experimental bounds on this mass 
exclude this possibility. We expect that a more general result, allowing for a large mixing 
angle, will not be so severely constrained, and enhancement of the cross section by as 
much as 50% (as in the left side of Fig. [3]) may be allowed. 

Analyzing the higher order corrections, in the MSSM we observe the same pattern as 
in the SM: large NLO corrections which further increase at NNLO. For rrisusy ^ 1 TeV 
the genuine supersymmetric corrections are suppressed and the difference from the SM is 
mainly due to the factor in Eq. flTHl) . However, the MSSM i^'-factors are very similar 
to their SM counterparts. This confirms the estimate of Ref. [2U] based on the complete 



NLO corrections to the matching coefficient and the NNLO calculation in the effective 
theory. 
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Our calculation has been performed in the limit of degenerate supersymmetric masses. In 
order to explore all interesting regions of the supersymmetric parameter space it would 
be necessary to relax this restriction. Since a three-loop calculation with many different 
mass scales is currently not feasible it is more promising to employ a strategy based on 
asymptotic expansions which has been applied successfully to predict the lightest MSSM 
Higgs boson mass in Ref. |26] . 
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